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1.0 


^5333SS3S3SSS=s* 

capabilities have been incorporated into IPOST to aid in PO^T-TlSSr 1616 ” and 

»“thTSal pl ““* and 

IPOST capabUities and 


feature 





CAPABILITY 


ptunization 
arameter* 


Maximum controls 
Control parameters* 



NPSOL 


AV magnitude, mass, time. 




nsitivity matrix 


Maximum events 
Event criteria* 


Event activities 


. -Ba 

Finite differencing, 
targeting 


00 




eed, c 


Maximum maneuvers/ 1 15 
subproblems I 


impul81ve AV « launch, orbit insertion, 


ettison 


Planetary bodies 


Trajectory 
rturbations 


> mi pu 

moon, any user-defined 




bodies* 116 planets ’ Earth ' 8 moon , any usei 

Analytic, precision (JPL) 

Central body, perturbing Dodies, radiation 
«JA aerodynapncs, thrust 




pressure, 


Table 1 - 1. IPOST Features/Capabilities 
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criteria, such as absolute a £^ g ’ activities can be initiated or 

consumption. At the event times, va , . ateerine law. changing 

terminated, such as employing a diffe performing an impulsive delta 

traiectorv propagators or propagation step size, penormmg uu v 

velocity maneuver or jettisoning a probe or stage. 

The time period between two ££S3to 

propagation takes place in each pba . dece a sin g computational speed): 

Oisted in order of l ^ rea81 ^ g E nXe Cowell Propagator selection depends upon 

tfSe .for parametric feasibility analysis, or 

precision detailed trajectories to support subsystem design. 

ssss 

optimization) for each search parameter value. 

SdSrn SlSdiSiSSt- uSJL ^parameters 

in^^nde^^arial^’o^contr^^mme ters^As ^rtof^or msteadof ^ ^ ^ 

function and the coX.fpaTa^emT 81 

values that meet specified target parameter conditions. 

Generalised targeting and optimisation uses ^ 

"s'that the trajecm^ <hrided into 

for gravity assist (swingby) missions. 
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easts an option to 

phase is divided into mdependenTXm»nL ^K . ^° d ' J° thi6 case ' each 

ag, ^»“ ssa: Sars** » 

time than explicit optimisation to ZcZve thS *"“**'' CP 

KffoV^ta^^^onSi' $TRA ^, and * TAB - *TOP contains a 

JTOP 1 Z D dlt S dat ® that descril >“ each lrission wnt^Cf ft 2 fu 
JTOP , and there must be one $TRAJ for each event STAR 14 ““P* follow 

tabular data such as thrust vs. time or drag c^°L! ^ B - S _ l ? ed t® mput 

angle of attack. Input and output units arem“ric “ aeh number and 
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2.0 ttFJ ATIONS WTP TO POST 

The software architecture of IPOST has been influenced by the* ^tecture of the 

F?Si? ing logic * general teble 

lookup capability, and the dictionary of output variables. 

The Input Processor accepts input in a namelist format from ^ input file and 

some Ihesenames to index pointers to the 

rt^Tv^s. SS.»« and processing of the data is also 
performed by the Input Processor. 

to POST). 

d^y^m^OST 6 events have been 

removed. 

Th*» evdina an d phasing have been restructured in accordance with trajectory 
^o2SoTloS? Sion phases have been grouped into sub-problems winch 
~ 3* sScific values Targeting may involve iterative computational 

^^u?h S S^"f-b-UlmpUe. 

each cycle to optimize the mission profile. Previously, all flight phases we 
computed. 

The design and structure for general table look up and interpoiation has b^n 
borrowedfrom the POST program. The dictionary of output variables has been 
borrowed from the POST program to support the design and capabilities described 
above Anv variable whose name appears in the dictionary may be written to the 
oS^ui dSs) used is the criterioSTor an event, a target in the decomposition 
sub^problem, a target of objective function of the mam program or as a control in 
either the sub-problem or the main problem. 
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3.0 PROGRAM MACROLOGIC Continued 



8 













4.0 SUBROUTINE HIERARCHY 


0 
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INTERPLANETARY POST CPC 


232 

Non-linear 
Programming - 
Stanford 
Optimization 
Laboratory 
NPSOL 
externals: 
CONFUN, OBJFUN 


1 Calculate values | 

and first 
partials of 

the 

non-linear 
constraints 
NPSOL; 
OQNFUN 

for 

2321 1 | 

1 Calculate values I 

and first 
partials of 
non-linear 

the 

constraints 

NPSOL; 

CNFWC 

for 


Calculate values 
and first 
partials of the 
objective 
function for 
NPSOL; 


2321 1 1 I 


Evaluate 

function 

(Perform 

Trajectory); 


TRAJM 


2321121 


Evaluate 
partlats of 
function 
analytically . 


ANPART 




2321 12 


Evaluate 

gradient of 

function 

(Perform 

Gradient of 

Trajectory); 

GRAD 


2321122 


Evaluate 
gradient of 
function wrt. 
independent 
variable (ks) 
(Perform 
Gradient of 
Trajectory); 
GRAD2 


232113 I 




2321 123 


Determine 
perturbation 
w/o loss of 
significance; 


Evaluate 
function 
(Perform 
T rajectory) ; 


TRAJM 



2321 1 23 


Determine 
perturbation 
w/o loss of 
significance 
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INTERPLANETARY POST CPC 


232111 


Evaluate 

function 

(Perform 

trajectory); 

TRAJM 

ALL MODULE NUMBERS BEGIN 2321 1 | 


..12 


Perform 

targeting 

this 

sub-problem; 

TRAJ 


..13 


Perform 
target this 
sub-problem; 


TTS 


..11 


Get elapsed 
CPU time 
in seconds; 


SECOND 


..14 


Terminate 
profil record 
and close 
profil file; 

CLSPFL 







INTERPLANETARY POST CPC 


2321112 







INTERPLANETARY POST CPC 


2321111 


ALL MOOULE NUMBERS BEGIN 232111 



^lUz 

..12 

..13 

L 

..14 

Tf 

Set input 

Iset initial 1 

Save initial 

Set independent 

Phase 

for the 

an 

conditions 

the values of 

initiali 

new 

for the new 

for the new 

independent 

ration; 

phase; 

phase; 

phase; 

variables into 
working storage; 

DINPT 

SETIC 

SAVIC 

SETIV 

PHZXM 


-16 


..1511 


J41J 

cr 

-151 | 

B i V 

1 

Cycling 

initiali- 

zation; 

USUBLK 

1 

| CVOOtl 



■ 

■■ 


1 

A 



..17 


Return target 
[and objective 
function values 
from the 
[trajectory; 

CALF 


.18 


Phase 

Terminate 

execu* 

profil 

tion; 

record t 


dose profil 


file; 

PHZXM 

CLSPFL 


-i5ia 


-1514 


Initialize 


Auxilary 


Initialize 


Initial 

ize 

trajectory 


caleula- 


fmeto 


output 

state and 


lions 


go: 


variables: 

propaga- 


initializa- 






tion; 


tion; 






OYNXMI 


AUXFMI 


TQOEM 


NFXMI 
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INTERPLANETARY POST CPC 
















4-0 SUBROUTI NE HIERARCHY Continue 

INTERPLANETARY POST CPC 


23211126 
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5.0 LIST OF IPO ST SUBROUTINES 


SUBROUTINE NAME PURPOSE 


ANPART 

ANPARTM 

ATMOS 

ATMOS2 

AUXFM 

BLKDT1 

CALF 

CDERIV 

CLSPFL 


Calculates the unweighted sensitivity 
matrix (CJAC) and the unweighted 
gradient of the objective variable (OBJGRD) 
using analytical methods supplied by the 
user. 

Calculates the analytic partials of the 
master problem dependent variables 
objective function with respect to the master 
problem independent variables. 

Calculates the pressure, density, 
temperature, and molecular weight at a 
given altitude for the primary planet. 
Currently, only earth data is available. 

Interpolates on the atmospheric tabular 
data to give the pressure, density, 
molecular weight, and temperature for a 
given altitude. 

Performs auxiliary calculations at the 
phase level. 

Block data initialization for targeting and 
optimization variables. 

Returns the value of the TARGET(I) and the 
objective function at event DEPPH(I) and 
OPTPH, respectively. 

Computes the derivative of the state with 
respect to time for use in integrating the 
equations of motion or the Cowell method. 

Terminates Profile record and closes Profile 
file. 


CNFUNC Evaluates nonlinear constraints and their 

derivatives v th respect to the variables 
U(I). 

COLLCI Initializes problem for the collocation 

algorithm 
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CONFUN 

COWELL 

CRAO 

CRADRX 

CYCXM 

CYCXMI 

DATA 


Computes the values and first derivatives of 
the nonlinear constraints. 

Performs the propagation of a trajectory 
using the Cowell method. 

Provides an indexing function to CRADRX 
to promote portability. 

Finds the addresses of the CRITR for each 
phase record, in IBKT, and finds the 

addresses of 

INDVR/DEPVR = 5HCRITR. 

Performs cycling functions and calculates 
the delta time step for the trajectory 
propagation. 

Performs cycling initialization. 

Block data initialization for real trajectory 
variables. 


DATA1 

DERIV 

DICT 

DICT1 

DIGDIF 

DENPT 

DTMDL 

DTMULT 

DYNSOl 


Block data initialization for integer 
trajectory variables. 

Computes derivative for the trajectory 
propagator. 

Block data initialization for names of real 
trajectory variables. 

Block data initialization for names of 
integer trajectory variables. 

Computes the number of significant figures 
in the difference A less B. 

Retrieves input data for this phase from 
GBKT and puts it in the trajectory variables. 

Modifies the value of DELT according to the 
user-selected simulation models. 

Calculates the variable step, DT, for the 
Multiconic propagation method. 

Trajectory propagator/dynamics 
integration for the Conic integrator. 
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DYNSl 


DYNS2 

DYNS3 

DYNS4 


Trajectory propagator/dynamics integration 
for the One-step integrator. 

Trajectory propagator/dynamics integration 
for the Multiconic integrator. 

Trajectory propagator/dynamics integration 
for the Cowell integrator. 

Trajectory propagator/dynamics integration 
for the Encke integrator. 


DYNXM 

Calls user-specified propagation method. 

DYNXMI 

Calls appropriate impulsive maneuver 
subroutine if applicable. 

EDERIV 

Evaluates the derivative of the deviation 
state matrix at any time. 

ENCKDT 

Determines the step size for the Encke 
propagator. 

ENCKE 1 

Performs trajectory propagation on the state 
of the S/C about the primary body with 
perturbations using the Encke method. 

ERRDVG 

Determines whether an error string 
diverges. 

ERROR 

Prints an error/waming message on the 
output device, then terminates the job if it is 
a fatal error. 

EXPN 

Input utility to create open space in the 
"bucket" by shifting existing data up into 
the array. 

FINDT 

Finds the normalized time value for 
Lambert's problem. 

FMATRX 

Evaluates the derivative of the acceleration 
applied to the S/C with respect to the 
position and velocity. 

GENTAB 

General table interpolation routine. 
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GRAD 

Calculates the gradients to each of the 
targets and to the optimization index with 
respect to the controls. 

GRAD2 

Calculates the gradient of each 
constraint/objective function with respect to 
the KSth control. 

GRADNR 

Calculates the sensitivity matrix with 
respect to perturbations of the controls. 

GRDCOL 

Performs the finite difference for subroutine 
GRADNR by the brute force method. 

HERM 

Determines the polynomial coefficients for 
the hermite cubic polynomial representing 
the state 

IMPULS 

Computes the impulsive delta velocity 
parameters for midcourse. 

INFXM 

Information output subroutine for 
controlling printing trajectory print blocks. 

INFXMI 

Performs quantization of the printed page 
print time and the profile print time. 

ISERCH 

Converts dictionary names to indices 
relative to IV. 

ISORTD 

Sorts an array in decreasing order using 
the Shell-Mezgar algorithm 

LAGRNG 

Fits a polynomial to three points using 
Lagrange's method. 

LAUNCH 

Calculates a hyperbolic escape trajectory 
from a park orbit including post-injection 
Cartesian state at periapsis. 

LOCF 

Provides an addressing algorithm for 
IPOST. LOCF is returned in units of 
INTEGER*4. 

MASSCK 

Checks for sufficient propellant mass for 
sub-events. 

MASSDT 

Mass integration for active thrust sub- 
events. 
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MASSET 

MASTER 

MD2050 

MFROMV 

MOTAB 

MOTION 

MPHI 

MPHIO 

MSRR1 

MSVPO 

MSVR2 

MSVRO 


General N-stage mass modelling 
ini tialization and jettison. 

Main program of the Interplanetary Post 
Simulation Program. 

Transforms a vector from mean equator 
nnH equinox 1950 (or mean 2000) to mean of 
date, or vice versa. 

Extracts a matrix from a given vector by 
columns. 

Returns the index to the greatest 
independent table-point less than or equal to 
the value of the independent variable. 

Trajectory state and propagation 
initialization. 

Computes the state transition matrix by 
chain rule from subproblem maneuver Ml 
to the next subproblem maneuver M2. 
(M2.GE.M1) and store in partitions. 

Computes the composite state transition 
matrix by the chain rule from one event to 
another. 

Computes the sensitivity matrix SRR.l, 
delta position (I) with respect to delta 
position (T-l) assuming intermediate target 
conditions are held constant for the pb 
maneuver. 

Computes the sensitivity matrix SVPO, 
maneuver velocity with respect to targets. 

Computes the sensitivity matrix SVR2, 
velocity at I with respect to position at 1-1 
assuming all other intermediate conditions 
are constant. 

Computes the sensitivity matrix SVRO, 
velocity at I with respect to position at I, 
assuming fixed target conditions. 


20 


MULTIC 


Performs the propagation of a trajectory 
from time TI by a delta time step of DT days 
using the STUMPFF- WEISS method. 


NCOUNT Determines the values of variables 

NSDEPV, NSINDV, and NSUB by counting 
INDXSD and INDXSI arrays. Also, 
determines the number of independent and 
dependent master level variables. 

NLOUT Outputs the summary of NPSOL's final 

solution. 

NLPRG Executive to implement NPSOL 

Optimization package into IPOST. 

NOMTAB Sets the table nodes to default values. 


NOMOUT 

NPRES 


Prints summary of NPSOL constraint/ 
objective and derivative evaluation. 

Computes constraint residuals. 


NPSLIC 

OBJFUN 

OFFBND 

ONESTP 

ORBINS 

PACC 

PAD 


Sets initial conditions for the Stanford 
University 

Optimization Laboratory Non-linear 
Programming Algorithm, NPSOL. 

Computes the value and first derivatives of 
the nonlinear objective function. 

Handles an out-of-bounds table exception. 

Performs the propagation of a trajectory 
from time TI to time TF using the One step 
method. 

Computes transfer maneuvers from an 
approach hyperbola to a desired elliptical 
orbit, including post injection state at 
periapsis. 

Computes perturbing accelerations for 
trajectory propagators. 

Determines the delta X that produces the 
most precise finite difference derivative 
without losing significance. 
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PAGER 

PARITO 

PBLOCK 

PHSTST 


Controls paging for output. 

Computes the master problem cost value 
and partial for the Interplanetary Targeting 
and Optimization Option (ITOO). 

Prints out various blocks of data as 
required. 

Performs test for abnormal phase 
termination. 


PHZCM 

PHZCMI 

PHZXM 

PHZXMI 

PLOT 


Performs the phasing functions for 
integration over a phase for collocation 

Performs phasing initialization 

Performs the phasing functions for 
integration over a phase. 

Performs phasing initialization. 

Collocation post processing 


PLOTIT 

FOLEQ50 


PRNPIC 

PRNTIP 

PSCAN 


Collocation post-processing 

Calculates the rotation matrix to convert 
from earth equatorial to planet equatorial 
system and vice versa. Also, calculates the 
transpose of this rotation matrix for 
rotating from planet equatorial to earth 
equatorial system. 

Prints initial user inputs to NPSOL 
Targeting and Optimization package. 

Prints the initial conditions for the 
trajectory and Targeting/Optimization 
algorithm for the current problem. 

Parameter scan master level target 
variables. 


READAT 


Input Processor Executive subroutine for 
the IPOST program. 
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RECTFY 

ROTRTN 

rotryp 

ROTUVW 

RTAB 

RTABX 

RTOP 

RTOPX 

RTRAJ 

RTRAJX 

RUNGE4 

SAVDAT 

SAVIC 

SAVSTT 


Determines if rectification is necessary, and 
if so, a new osculating orbit is found by 
choosing a starting point that coincides 
with the true orbital path. 

Crates the rotation matrix from the 
inertial equatorial system to the RTN 
system and vice versa. 

vp^o^ at€ fl th il? >tati0n matrix to convert a 
J ,?° r J? system to a vector in 

P™ 8y8tem > the UVW system, 
ttie merfa^ equatorial system, or the 
inertia I ecliptic system. 

Calcuiates the rotation matrix from the 
inertial equatorial system to the UVW 
system and vice versa. 

Reads Namelist $TAB. 

K;r brouttoe for reading the taMe 

feads targeting and optimization Namelist 

Executive subroutine for reading targeting 
“ PU ‘ and St0ring H “ ‘he 

Reads trajectory Namelist input. 

Executive subroutine for reading trajectory 
input and storing it in the data "bucket". 

Performs fourth order Runge-Kutta 
integration on the state. 

Saves default AAIVC to ZZENDC (Block 
data) values over multi-problem execution. 

the current state of the trajectory 
(event -) for quick function evaluation * 

£3Ss ne/ optimizati ° n 

Saves the state. 
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SCAREA 

SEARC 

SECOND 

SERCH 

SETCNT 

SETIC 

SETIV 


SETSTT 

SETWT 

SNMOUT 

SOICHK 


SRA02 

SRADP 

SRADRX 

TARG1S 


Calculates the S/C area normal to the 
velocity vector through the atmosphere. 

a q! rmrlfVlTTlR 


Gets the elapsed CPU time in seconds. 

Converts dictionary names to indices. 

Counts the number of independent and 
dependent variables in each sub-problem 

Restores the trajectory state previously 
saved by SAVIC. 


Sets the Targeting/optimization 
independent (control) variables to the 
desired values, based on the calculated 
control correction. 


Resets the state. 

Sets weighting for Newton-Raphson 
targeting method for sub-problems. 


Writes out subproblem status. 


Verifies if given radius vector in primary 
body reference frame is still within the 
primary body’s sphere of influence. If it is 
not, change state vector and primary body to 
reflect which sphere of influence the radiu 


Provide an indexing function to CRADEX to 
promote portability. 

Computes solar radiation pressure for 
trajectory propagation. 

Processes the namelist $TOP variables. 
Reorders the independent and dependent 
variables. 


Special One step targeting option which 
varies pseudostate to achieve desired target 
conditions. 
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TARGET 

TBADRX 

TGOEM 

TGOEMI 

TPRP 

TRAJ 

TRAJC 

TRAJM 

TSPXM 


Converts state vector into specified conic 
parameters and impulsive delta velocity. 
Sets addresses for tables; converts names of 
table arguments to addresses. 

Computes time-to-go to the next pending 
event. 

Time-to-go initialization. 

Finds the closest approach and time of 
closest approach of the vehicle to a 
secondary body for a given time interval and 
a primary central body. Only the attraction 
of the primary body is considered. 

Propagates the trajectory phase-to-phase 
over one subproblem. 

Performs logic to evaluate the collocation 
trajectory from phase, iensi, to phase, iensf, 
phase-to-phase. 

Performs logic to propagate trajectory 
subproblem to subproblem. 

Top level simulation module. 


TTS Target this subproblem. Propagates a 

trajectory to achieve a desired targeting. 

USUBLK Sets the USUB local jacobian matrix for sub- 

problem control. 
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6.0 LIST OF UTILITY SUBROUTINES 


E-ENTRY POINT 
F-FUNCTION 

NAME S-SUBROUTINE 

VARIABLE 

ANBET (A,B) 

F 

Function ANBET 
calculates the smallest 
angle between any two 
vectors A and B. 

ANOMLY (E,IOPT,AM,F) 

S 

Subroutine ANOMLY 
converts true anomaly 
to mean anomaly or 
vice versa. 

CAUUL (ICODE,DAY,IYR, 

MO, IDAY,IHR,MIN,SEC) 

S 

Subroutine CAUUL 
converts Julian date to 
calendar date or 
vice versa. 

CARCON (GMU, STATE, ORBT) 

E 

Entry Point CARCON in 
Subroutine CHGELM 
calculates orbital 
elements from a state 
vector. 

CHGREF (T,NOLD,NNEW, 
STATEO, STATEN) 

S 

Subroutine CHGREF 
converts a state vector 
from one celestial 
body centered 
ecliptic reference frame 
to a new celestial body 
centered ecliptic 
reference frame. 

CONCAR (GMU, STATE, ORBT) 

E 

Entry Point CONCAR in 
Subroutine CHGELM 
calculates the state 
vector from the orbital 
elements. 

DANBET (A33EF) 

F 

Function DANBET 
calculates the directed 
angle between any two 
vectors A and B. 
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NAME 


E-ENTRY POINT 

F-FUNCTION 

S-SUBROUTTNE 


VARIABLE 


EQTECT (DJ,XIN,XOUT) S 

ECTEQT (DC.XIN.XOUT) E 


EQ5EC5 (DJ,XIN,XOUT) 


EC5EQ5 (DJ,XIN ( XOUT) 


EPH50 (N,T,STATE) 


EPHEMS (NP,T,IEPHEM, STATE) S 


EPHEMP (T,IBODYl,IBODY2, 
PSTATE) 


Subroutine EQTECT 
transforms a vector from 
equator of date to ecliptic 
of date. 

Entry point ECTEQT in 
Subroutine EQTECT 
transforms a vector from 
ecliptic of date to equator 
of date. 

Entry Point EQ5EC5 in 
Subroutine EQTECT 
transforms from equator 
of 1950.0 to ecliptic of 
1950.0. 

Entry Point EC5EQ5 in 
Subroutine EQTECT 
transforms from ecliptic 
of 1950.0 to equator of 
1950.0. 

Returns planet position 
and velocity vectors (R,V) 
in mean ecliptic and 
equinox of 1950.0. 

Subroutine EPHEMS 
finds the state vector of a 
specified body at a 
specified time in 
heliocentric mean 
ecliptic and equinox of 
1950.0. 

Calculates position and 
velocity of an 
EPHEMERIS body from 
JPL Chebyshev 
coefficients. 
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NAME 


E-ENTRY POINT 
F-FUNCTION 
fi.ST TBROUTCNE 


VARIABLE 


EPHSAT (T,XPRI,N ,IPRIB , 
ISATB.XSAT) 


GOODYR (XI.TAU.GMU, 
IPFLAG,PSI,XF,PHI, 
PHINV) 

IFRM50 (SIN, SOUT) 


JACOBI ( JTARG.X, AUX.Y ,S) S 


Subroutine EPHSAT 
finds the state vector of a 
satellite orbiting its 
primary body at a 
specified time. 

Subroutine GOODYR 
solves a two body problem 
with partial derivatives. 

Subroutine IFRM50 
converts the input state 
vector from the I/O body 
reference frame to 
ecliptic mean of date 1950. 

Subroutine JACOBI 
computes 

the 3x6 sensitivity matrix 
of S = dY/dX by finite 
differencing. JACOBI 
presumes the subroutine 
TARGET has been called 
to compute and store the 
rt nminal Y parameter. 
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name 


E-ENTRY POINT 

F-FUNCTION 

S-SUBROUTINE 


VARIART.E 


LAMBRT (R1,R2,TFS,NREV, 
GMU.V1.V2) 

S 

Subroutine LAMBRT 
solves LAMBERTS 
problem by 

determining an orbit 

given two position 
vectors and a interval. 

LUBKSB (A,N,NP,INDX,B) 

S 

Subroutine LUBKSB 
solves the set of N linear 
equations AX = B. Not as 
the matrix A but rather 
as its LU decomposition 
determined by LUDCMP. 

LUDCMP (A,N,NP,INDX,D) 

S 

Subroutine LUDCMP 
replaces 

an NxN matrix A with 
the lower-upper (LU) 
decomposition of a row 
wise permutation of 
itself. 

M50IFR (SIN.SOUT) 

S 

Subroutine M50IFR 
converts an input state 
from ecliptic mean of 
data 1950 to I/O body 
reference frame. 

MADD (A,SA,B,SB,N,M,C) 

S 

Subroutine MADD 
performs matrix addition 
and scalar 
multiplication. 

MATOUT (A,N,M) 

S 

Subroutine MATOUT 
prints an NxM matrix A 
to output device. 

MCOPY (A,N,M,B) 

S 

Subroutine MCOPY 
copies matrix A to 
matrix B. 

MIDENT (NA) 

S 

Subroutine MIDENT 
transforms the matrix A 
into an identity matrix. 
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NAME 

E-ENTRY POINT 

F-FUNCTION 

S-SUBROUTINE 

VARIABLE 

MINV (A,N,B,ISING) 


S 

Subroutine MINV 
calculates the inverse, if 
it exists, of matrix A. 

MMULT (A,B,L,M,N,C) 


S 

Subroutine MMULT 
performs 

matrix multiplication. 

RDCHBY (ET.IBODY.CHBYFL, 
NCJPARAM,ICENT) 

S 

Reads Ephemeris file 
created by PREPRO for 
the current body and 
time. 

RELATV (T,NPRIM,NSEC,RHO) 

S 

Subroutine RELATV 
finds the state vector 
representing the position 
and velocity of a 
secondary celestial body 
with respect to a 
specified central celestial 
body. 

TARGET (ITARG,X,GMU,Y, 
AUXO) 


S 

Subroutine TARGET 
converts a state vector 
into specified conic 
parameters. 

VANGO (GAM,XNODE,XINC, 
VECTOR) 

S 

Subroutine VANGO 
computes a unit pointing 
vector. 

VCOPY (A,N3) 


S 

Subroutine VCOPY 
copies vector A to 
vector B. 

VCROSS (A3,C) 


S 

Subroutine VCROSS 
computes the cross 
product of two vectors. 

VDOT (A,B) 


F 

Subroutine VDOT 
computes the dot product 
of two vectors. 
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NAME 

E-ENTRY POINT 

F-FUNCTION 

S-SUBROUTINE 

VARIABLE 

VMAG(A) 

F 

Subroutine VMAG 
calculates the 
magnitude of vector A. 

VSCALE (S,A,N,B) 

S 

Subroutine VSCALE 
scales a vector A and 
outputs vector B. 

VUNIT (A3) 

S 

Subroutine VUNIT 
calculates unit vector B 
from vector A. 

VZERO (N,A) 

S 

Subroutine VZERO 
converts vector A to a 
zero vector. 
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7.0 GLOBAL COMMON BLOCKS 


COMMON AAIBKT - 

COMMON AAIVC - 
COMMON ATMOSC - 

COMMON CYCNC - 

COMMON CYCVC - 

COMMON DVIVC - 
COMMON DYNVC - 

COMMON EPHEJC - 

COMMON ER 
COMMON EPHEMC - 

COMMON GENIC - 
COMMON GENRL - 

COMMON INFVC - 
COMMON INFV1C - 

COMMON INFVC - 
COMMON IPRTO 

COMMON ITOOJC - 
COMMON MNVRJC - 


Data structure containing event definition records, 
tables, and phase state transition matrices. 

Defines beginning of the dictionary. 

Contains real type variables defining the atmospheric 
data for the earth. 

Con tains integer type variables used to perform cycling 
functions. 

Contains real type variables used to perform cycling 
functions. 

Contains delta velocity variables. 

Contains time references and other variables used to 
perform dynamics functions. 

Contains character type variables defining the planets’ 
alphanumeric names. 

Internal Table errors common. 

Contains the variables for the analytic planetary 
ephemeris. 

Contains miscellaneous top level variables. 

Data structure containing real and integer type input 
from namelist $TRAJ. 

Contains printout control variables. 

Contains the title and profile file identification 
character variables. 

Contains miscellaneous top level input control variables. 

Contains internal variables used to control printout 
during targeting and optimization. 

ITOO input variables common. 

Contains integer and character type variables for 
impulsive maneuver sub-events. 
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COMMON MOTBLC 
COMMON PACVC 

COMMON PADC 

COMMON PADSWC 
COMMON PAGERC 
COMMON PHZNC 

COMMON PHZVC 

COMMON PLANEC 
COMMON PROPDC - 

COMMON PROP JC - 

COMMON REDAT - 
COMMON SEARC - 

COMMON SENSIT - 

COMMON SERVC 

COMMON STATEC - 

COMMON STATNC - 

COMMON SUBPVC - 

COMMON TGONC - 


in 


- Table Nodes common (also MOTBLO). 

- Contains perturbing accelerations for trajectory 

propagation. J 

- Contains variables used in controlling perturbation 

tne targeting and optimization algorithm. 

- Contains internal sub-problem control variables. 
Contains variables used in pagination control. 

WwV 11 ^ 61 ^ t31>e variables associated with phasing 

^tionV 681 ^ Variables associated with phasing 

Contains planet related variables. 

Contains real variables controlling trajectory 
propagation and perturbing influences. 

Contains integer variables controlling trajectory 
propagation and perturbing influences. 

Contains input processing variables. 

vari^S ma i ter 1 j Vel feting and optimization 
variables, independent and dependent variables. 

T trix ^ration parameters for 
swingby targeting and op timiz ation. 

arrays 6 C0mm0n containin & constants and temporary 

orWwSemenfe' 16 ' 8 P ° Si ‘ i0n “ d Velodty and *• 

ch ? ra< * er type variables defining 
tne vehicle s position and velocity. ^ 

Contains sub-problem level targeting variables 
independent and dependent variables. 

Contains integer type time-to-go variables. 
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COMMON TGOVC 
COMMON ZZENDC 


Contains real type time-to-go variables. 
Defines the end of the dictionary. 
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8.0 



block data blkdtl 

***Wj dtreaI *8 (a-h,o-z) 
c*** blkdat 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


' viauat 

-*** program values 

the data statements in this routine are standardized 


standard fortran declarations 
common 
dimension 
equivalence 

type (can be anywhere above) 
data (can be anywhere) 
namelist 


include 'aaibkUnc' 
include 'genic.inc' 
include inflame' 
include 'inpvc.inc' 
include 'itoojc.inc' 
include padc.inc' 
include ’pagerc.inc’ 
include ’parmsc.inc' 
include 'prOfc.inc' 
include 'searc.inc' 
include ’servc.inc' 
include ’8ubpvc.inc' 


parameter (mxinv5=5*imrinv) 
c parameter (°«>nv2=3‘manVm>dnv + 41*nm„ v ) 


c — . 

C include 'aaibkt.inc' 



c common 

data 

1 ihkt(iyo 

c 

/aaibkt/ 

/ 






c include 'genic’ 

. . 

C common 

data 

1 iprt /63 
c 

: — — 

/genic / 

/ 



as 



C include ’infjc’ 


common /infjc / 
data 

1 iprofl / 0 / 

data 

1 npvar / 33 / 


c 

c 


common /infjcc / 

data outnam 

0 /'time ' , , 

1 ’x ”v \'z *, vx ,vy ,VZ 

2 ’sma ’ y ','eccen ','inc VanlongVargp , meaan 

3 ,'truan ','rperi ','vperi 

4 i'radiua', ’speed Vbtheti vinfy. , vmfa 

5 .’bdti ','bdri Vhypta ,c3i , ra i, ’ d ‘ hdr0 • 

6 .'btheto',’ v-infxo, vmfyo vinfzo , bdto .bdro 

7 ,’hypta ’,’c3o , Vrao Vdeco Jpa , alft 

8 .longp ’,’lat .long ,1) » 53 


include 'inpvc’ 



common 

data 

/inpvc / 

1 

ioflag 

/3 / 


common 

data 

/inpcc / 

1 

iprops 

data 

/maxevt*'none’/ 

1 

mnvrs 

/maxevt*’none7 


include ’itoojc’ 


common 

data 

1 itoowt 


/itoojc / 
/ ’sumass’ / 


c 

c 

c 

c 

c 


1 

2 

3 


include 'padc' 
common /padc / 


d npad / 9-dO, 4.d0, 14.44943979d0 

,pdlmax / 2.0d0 / 

.sigdel / O.OdO / 


/ 
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4 ,sdif / O.OdO / 
c 


c 

c 

c 


include ' 

pagerc 

c 


common 

data 

/pagerc / 


1 

icase 

/O / 


2 

,nlines 

/60 / 


3 

.header 

/" / 


4 

.npage 

/I / 

c 


common 

data 

/pager2 / 

c 

c - 

1 

ptime 

/" / 

c 

c - 


include 'prOfc.inc' 

c 





data idxdat / prsize * 0.0 / 
datahdrdat /prsize*' 7 
data frsttl / 0/ 
data nxtrec / 3/ 
data dm / 1/ 
data cdm / 2/ 
data chm / 2/ 
data dep / 34/ 
data cdep / 0/ 
data chep / -1/ 
data tdpos /0/ 
c 

c — 
c 

c — 
c 

1 

2 

3 

4 

6 

7 

8 

9 

0 

1 

2 

3 

4 


include 'searc.inc' 


common /seric / 
data 

ideb /0 / 

,ifdeg /25 * 0/ 

!2^22, 23,24, 

gi9^>a 24 2 / ^ 2 ’ 3 ’ 4,5 ’ 6,7 ’ 8A10>11 ' 12 > 13 ' 14 - 15 . 16 .lV,18,19,20 

,ipro /0 / 

.isrndep /0 / 

,ismind / lmxinv*0 / 

»mxitop /10, 5*-l/ 

,ndepv /0 / 

,nindv /0 / 

,nclin /O / 


37 



5 .nsminc /O / 

data 

1 nsmind /O / 

2 ,opt /O / 

3 ,pert /nurinv * .OOOldO / 

4 .smincr /l.dO / 

5 ,8rchm /'none' / 

6 .isens /0 / 

7 ,etanl /.9d0 / 

8 ,pmtpd /.false./ 

9 ,sclflg /.false./ 

0 .wgtflg /.true. / 

1 .ftol /l.d-14 / 


2 

indplb /mxinv*-10.dl0 / 



3 

,indpub 

/mxinv^+lO.dlO 

/ 


data 




1 

depvlb /mxinv*-10.dl0 / 



2 

.depvub 

/mxinv*+10.dl0 

/ 


3 

.deptl 

/mxinv*1.0d0 

/ 


4 

.iter /0 

/ 



5 

.nstep 

/0 / 



6 

.inform 

/0 / 



7 

.nrowa 

/10 / 



8 

.nrowj 

/7000 / 



9 

.nrowr 

/7000 / 



0 

.nctotl 

/0 / 



1 

.msglvl 

/10 / 



2 

.nstate 

/0 / 




data 




1 

istate 

/mxinv*0, mxinv*0 

/ 

2 

.clamda 

/mxinv*0.dO, nmnv*0.d0/ 

3 

.orthog /.true. / 



4 

.leniw 

/mxinv5 / 



5 

,lenw 

/mxinv2 / 



6 

.icnfun 

/0 / 



7 

.nbase 

/2 / 



8 

.ndigit 

/56 / 



9 

.epsmch 

/2.78d-17 / 



0 

.flmin 

/2.94d-39 / 



1 

.flmax 

/1.70d+38 / 



2 

,nin /5 

/ 



3 

.nout /6 

/ 



4 

.rteps 

/l.d-7 / 



5 

.epsaf 

/l.d-14 / 




6 .weight /mxinv*0y 


c 

c common /sercc / 
data 

1 istm / 'forward' / 


38 


2 npinpt / npinput' / 


common /serfc / 
data 

1 fesn /999.0d0 / 


c common /tofioc / 

c data de P vl /mxinv * O.OdO/, optvl/0.0d0/,u/mxinv*0.0d0/ 

c common /tofeic / 

data 

1 depph /mxinv * 900.0d0 / 

2 ,depvr /mxinv * O.OdO / 

3 ,indph /mxinv * 0 / 

4 ,indvr /mxinv * 0 / 

5 .optph /900.0d0 / 

6 .optvar /O.OdO / 

7 .sindph /mxinv*0.d0 / 

8 .namsvr /mxinv*0 > mxinv*0/ 

9 ,namopt / 0/ 

0 ,wvu / mxinv*0.d0 / 

1 ,wvlc /10*1.0d0 / 

2 ,wvnlc /mxinv*0.0d0 / 

3 ,wopt /l.OdO / 


include 'serve' 


common /serve / 

data 

L temp /50*0.d0 / 

,stemp /25*0.0d0 / 

1 *Pi°2 /L57079632679489d0 
• .Pi /3.1415926535898d0 
! '^d /0.0174532925199433d0 
.dpr /57^95779513082d0 
,twopi /6.2831853071796d0 
.spdy / 86400.d0 / 

,sphr / 3600.d0 / 

,fpl /l.OdO / 

,fpem6 / l.0d-6 / 

,fpem8 / i.0d-8 / 

,fpeml0 / l.0d-10 / 

,fpeml2 / 1.0d-12 / 

,idebug /O / 

,iephem /0 / 

data 

unitm /l.dO ,0.d0 ,0.d0 
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CO M M 


,0.d0 ,l.dO ,0.d0 
,0.d0 ,0.d0 ,1-dO 
/ 

data 

zerom /O.dO ,0.d0 ,0.d0 
,0.d0 ,0.d0 ,0.d0 
,0.d0 ,0.d0 ,0.d0 
/ 

data 

1 enOis / l.Od-S / 

2 ,xinf / 1.0d+38 


c 

c 

c 

c 

c 


include subpvc' 


common 

data 


/subvpc 


/ 


1 

2 

3 

4 

5 

6 

7 

8 
9 

0 

1 
2 

3 

4 

5 


atarwt / 2.d0, 7.d0, 3.d0, l.dO, l.dl, 8.d0, 
5.d0, 0.7d0, 2*0.d0 / 

,indslb /mxnsvr*-l.d38/ 

4n dsub /mxnsvr*l.d38/ 

,indxsd / mxnsvr*0 / 

,indxsi / mxnsvr*0 / 

,isub / mxnsub*0 / 


,nudtar 
,npi /7 

,nsdepv 

^sindv 

,nsub 

,pertsb 

,targmt 


/I / 

/ 

/0 / 

/0 / 

/o / 

/ mxnsvr*0.0001d0 
/ mxnsub*1.0d0 / 


„ ,rfsoi / mxnsub*0.5d0 

6 ,optsvT / mxnsub*0.d0 / 

7 ,optsph / mxnsub*900 / 

8 ,tolf / l.dO/ 

9 ,tolu / l.dO/ 

common 
data 

1 modelt 
end 


/ 


/subpcc / 

/ mxnsub*’nraph 
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block data data 
implicit real*8 (a-h,o-z) 

g************************************************************* 

c 

c data - defines computational commons 
c (aaiv - zzend) computational values 
c 

c the data statements in this routine are standardized 
c 

g************************************************************* 

c 

include 'aaivc.inc' 
include 'atmosc.inc' 
include 'cycvc.inc' 
include 'dvivc.inc' 
include 'ephemc.inc' 
include 'infvc.inc' 
include ’motblc.inc' 
include 'pacvc.inc' 
include ’phzvc.inc' 
include 'planec.inc' 
include ’propdc.inc’ 
include 'sensit.inc 
include 'statec.inc' 
include 'tgovc.inc' 

c 

c *********************** ************* *4,4, ********************** 

c 

c *** computational data 

c 

g***********************************,^,^^^^^^^^^^^ 

c 

c include 'aaivc.inc' 


c common /aaivc / 

data 

1 iv(0) /0 / 

c 

c *** common atmosc 
c 

data ali/0.d0,11.019d0,20.063d0,32.162d0,47.350d0,52.43d0,61.59d0, 
+ 80.d0,90.d0,100.d0,110.d0,120.d0,150.d0,160.d0,170.d0, 

+ 190.d0,230.d0,300.d0,400.d0,500.d0,600.d0,700.d0, 

+ 28*0.d0/ 

dataatms/1.225dO,1.01325d5,28.964dO,3.139d-7,2.87d2,7.12dO, 

+ 0.d0,1.4d0,0.d0,0.d0/ 

data tpt/288.1d0,216.65d0 > 216.25d0 > 228.65d0,270.65d0,270.65d0, 

+ 252.65d0,180.65d0,180.65d0^10.65d0,260.65d0360.65d0, 
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+ 960.65d0,1110.65d0,1210.66d0,1350.65i0,1550.65d0, 

+ 1830.65d0,2160.65d0,2420.65d0,2590.65d0,2700.65d0, 

+ 28*0.d0/ 

An ts* pres/l.d0,2.284d-l,5.462d-2,8.567d-3,l-095d-3,5.823d-4, 

+ 1.797d-4,1.024d-5,1.622d-6,2.98d-7 ,7 .22d-8,2.488d-8, 

+ 5.0d-9,3.64d-9,2.756d-9,1.66d-9,6.869d-10,1.86d-10, 

+ 3.977d-ll,108d-ll,3.4d-12,1.176d-12 t 28*0.dCV 

data wgi/9*28.964dO,28.88dO,28.56dO,28.07dO,26.92dO T 26.66dO, 
+ 26.5d0,25.85d0,24.69d0,22.66d0,19.94d0,17.94d0,16.84d0, 

+ 16.17d0,28*0.d0/ 


*** 


cycvc common 


data date / O.dO, O.dO, O.dO, O.dO, O.dO, O.dO/ 

data dt / 1.0d+6/ 

data dtime / O.dO / 

data dtimr/ l.dO, 19 * O.dO / 

data dtO / O.dO / 

data juldat / O.dO/ 

data tdurp / O.dO/ 

data time / O.dO/ 

data timrb / 20 * O.dO/ 

data timrf / 20 * O.dO/ 

data tref / O.dO/ 

data trefp / O.dO/ 

c 

c *** dvivc common 


data dvx / 3*0.d0/ 
data dvmag / O.dO/ 
data dvsum / O.dO/ 


c *** ephemc common 
c 

data gmu / 1.327124419330078d+ll 
m ^ .20320 8080433758d+4 

v ,3-248587705064894d+5 

e , 3.986004 850429589d+5 

m ,4-28^28685339965d44 

j ,1.266858077864550d+8 
s ,3.792449009286876d+7 

u , 5.8020 58945865924d+6 

n ,6.850305659666957d+6 

p ,4.423748064433594d44 

m ,4.9027927809d+3,12*0.d0 

z/ 


c 

c 

c 


radius of the sphere of influence 
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data reoi / 1.000000000000d+ll 


m 

,3. 18987 802284 ld+05 

V 

,1.45833 656623 3d+06 


e 

,2.167226195872d+06 


a 

,1.563493250956d+06 


j 

,7.66407843 1 145d+07 


s 

,9.399359395396d+07 


u 

,1.0100 8450 0916d+08 


n 

,1.677463630809d+08 


P 

,5.000 0000 00000d+08 


m 

,6.6d+04,12*0.d0/ 


time radius of sphere of influence 

data tsoi / 1.00d+ll 


m 

,2.00d+01 


V 

,2.00d+01 


e 

,2.00d+01 


a 

,2.00d+01 


j 

,7.00d+02 


6 

,7.00d+02 


U 

,7.00d+02 


n 

,7.00d+02 


P 

,2.00d+02 


m 

,0.5d0,12*0.d0/ 



semi-major axis time polynomial coefficients 


of solar-system bodies 
data 


z 

((cGj, 1)4=1,4)4=0,22) 


s 

/4*0.d0 


m 

,5.790977169d+07 

,3*0.d0 

V 

, 1.0820988790d+08 

,3*0.d0 

e 

,1.495998116d+08 

,3*0.d0 

a 

,2.279418963d+08 

,3*0.d0 

j 

,7.783076165d+08 

,3*0.d0 

8 

,1.429411131d+09 

,3*0.d0 

U 

,2.875072804d+09 

,-5.568 

$ 

,2*0.d0 


n 

,4.504503336d+09 

,-24.89 

$ 

,2*0.d0 

P 

,5.890213800d+09 

,3*0.d0 

m 

,3.843984402d+05 

,3*0.d0 

$ 

,48*0.d0 


z 

/ 



eccentricity time polynomial coefficients 

of solar-system bodies 



data 


z 

((c(ij, 2),i=l, 4)j=0, 5) 
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6 

/4*0.d0 


m 

,.20563175 ,2.0406d-5 

$ 

,-2.84d-8 

,-1.7d-10 

V 

,6.77188d-3 

,-4.7766d-5 

$ 

,9.75d-8 

,4.4d-10 

e 

,.01670862 

,-4.2037d-5 

$ 

,-1.236d-7 

,4.d-ll 

a 

,.09340062 

,9.0483d-5 

$ 

,-8.06d-8 

,-3.5d-10 

j 

,.04849485 

,1.63244d-4 

$ 

,-4.719d-7 

,-1.97d-9 

z 

/ 



data 


z 

((c(ij,2),i=l, 4)j=6, 22) 


8 

/.05550862 

,-3.46818d-4 

$ 

,-6.456d-7 

,3.38d-9 

U 

,.04629590 

,-2.7337d-5 

$ 

,7.9d-8 

,2.5d-10 

n 

,8.98809d-3 

,6.408d-6 

$ 

,-8.d-10 

,-5.d-ll 

P 

,2.488033d-01 

,3*0.d0 

m 

,5.4900489d-02 

,3*0.d0 

z 

,48*0.d0 


z 

/ 



inclination time polynomial coefficient] 


of solar-system bodies 
data 


z 

((c(ij,3),i=l,4)j=0,5) 


8 

/4*0.d0 


m 

,.122260070 

,-1.0387d-4 

$ 

,1.414d-8 

,7.156d-10 

V 

,.05924803 

,-1.4954d-5 

$ 

,-5.662d-7 

,1.75d-10 

d 

,0.d0 

,-4.2078: 

$ 

,7.0983d-7 

,-2.316d-8 

a 

,.0322838 

,-1.42208d-4 

$ 

,-3.936d-7 

,-4.71d-10 

j 

,.02274635 

,-3.4683d-5 

$ 

,5.791d-7 

,1.6d-9 

z 

/ 


data 


z 

((c(io,3),i=l,4)j=6,22) 


6 

/. 043439 12 

,4.4532 ld-5 

$ 

,-9.557d-7 

,3.14d-10 

U 

,.0134948 

,-2.94496d-5 

$ 

,6.0912d-8 

2-8d-10 

n 

,.03089149 

,3.9392d-6 

$ 

,4.014d-9 

,0.d0 


£996707636393978d-01 


,3*0.d0 
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m 

,8.9804108d-02 

,3*0.d0 

z 

,48*0.d0 


z 

/ 



longitude of ascending node polynomial coefficients 

of solar-system bodies 

data 

z 

((c(ij,4),i=l,4)j=0,5) 


s 

/4*0.d0 


m 

,.84353321 

,-2.189026d-3 

$ 

,-1.5416d-6 

,-3.421d-9 

V 

,1.3383171 

,-4.852155d-3 

$ 

,-2.4881d-6 

,-3.456d-9 

d 

,3.0521127 

,-4.207828d-3 

$ 

,7.0983d-7 

,-2.316d-8 

a 

,.86495189 

,-5.44845d-3 

$ 

,-1.11689d-5 

,-3.7402d-8 

j 

,1.753435277 

,3.083697d-3 

$ 

,1.57755d-5 

/ 

data 

,-1.277d-7 

z 


z 

((c(ij,4),i=l,4)j=6,22) 


8 

/1.98383764 

,-4.47965d-3 

$ 

,-3.2018d-6 

,6.23d-9 

U 

,1.291647441 

,1.29409d-3 

$ 

,7.07556d-6 

,1.82d-9 

n 

,2.300065696 

,-1.07601d-4 

$ 

,-3.82d-8 

,-1.36d-9 

p 

,1.914337511979504d+00 ,3*0.d0 

m 

,4.52360 1515d+00 ,-9.2422d-04 

$ 

,3.6267d-05 

,3.4d-08 

z 

,48*0.d0 


z 

/ 



argument of perigee polynomial coefficients 

of solar-system bodies 

data 

z 

((c(ij,5),i=l,4)j=0,5) 


8 

/4*0.d0 


m 

,.50833109 

,4.961748d-3 

$ 

,1.307d-6 

,4.102d-9 

V 

,.957902789 

,4.93706d-3 

$ 

,-2.16379d-5 

,-8.961d-8 

d 

,-1.255517014 

,9.83749d-3 

$ 

, 1.9127 ld-6 

,3.1503d-8 

a 

,5.087672137 

,.0128958 

$ 

,8.1458d-6 

,4.2638d-8 

j 

,-1.503306748 

,6.784042d-4 

$ 

,-3.165154d-6 

,4.2621d-8 

z 

/ 



data 
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8 

$ 

U 

$ 

n 

$ 

P 

m 

$ 


((c(ij,5),i=l,4)j=6,22) 
/- .3596903 
,1.2419d-5 
,1.7278622 
,-8.728d-6 
,-1.46014884 
,1.2688d-6 


,.0143678 
,7.898d-8 
3-64845d-4 
,5.39d-9 
,6.1652d-4 
,9.6d-10 


,1.995581715995208d+00 ,3*0.d0 
,1.311550025d+00 3.868587d-03 
,-2.16472d-04 ,-2.43d-( 


z 

,48*0.d0 


z 

/ 



mean anomaly polynomial coefficiec 

of solar-svstem bodies 

data 

z 

‘ (ij.f l,4)j=0,5) 


8 

/4* 10 


m 

3.0 07445 

3608.7875 

$ 

,4.5691d-3 

,-6.46d-10 

V 

,.879926836 

,1021.32847 

$ 

,2.415483d-5 

,9.303d-8 

d 

,-.043125323 

,628.3019552 

$ 

,-2.72167d-6 

,-8.34d-9 

a 

,.3381233516 

,334.0534958 

$ 

,3.06864d-6 

,-5.288d-9 

j 

,.3494179706 

,52.9653341 

$ 

,-1.4094d-5 

,8.018d-8 

z 

/ 

data 

z 

((cGJ,6),i=l,4)j=6,22) 


8 

/-.750 1305745 

31.32002138 

$ 

,-5.55102d-6 

,-8.554d-8 

U 

3.461784222 

,7.476600918 

$ 

,1.568d-6 

,-7.103d-9 

n 

,4.471969422 

,3.812794648 

$ 

,-1.22d-6 

3.7d-10 


P 

$ 

m 

$ 

z 

z 


,2*0.d0 

,4.719966573d+00 

,-1.9774d-05 

,48*0.d0 


3-29971481d-01 

,3.3d-08 


c vector right ascension polynomial coefficients 

c of solar-system bodies 
data 

z ((polev(ij,l),i=l,4)j=0,22) 
s /286.0193d0 ,3*0.d0 

m 380.86554d0 ,-.03289d0 
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$ 

,.0000 ldO 

,0.d0 

V 

,98.02255d0 

,3*0.d0 

e 

,-1.3435d-6 

,-.6402780091d0 

$ 

,-8.39481d-5 

,-5.0003d-6 

a 

,316.8538d0 

,-.0996d0 

$ 

,2*0.d0 


j 

,268.0447d0 

,3*0.d0 

8 

,38.41314d0 

,3*0.d0 

U 

,76.761d0 

,3*0.d0 

n 

,295.57 12d0 

,3*0.d0 

P 

,313.89136d0 

,3*0.d0 

z 

,52*0.d0 


z 

/ 



pole vector declination polynomial coefficients 

of solar-system bodies 

data 

z 

((polev(ij,2),i=l,4)j=0,22) 

8 

/63.7718d0 

,3*0.d0 

m 

,61.39767d0 

,-.0047 ldO 

$ 

,.0000 ldO 

,0.d0 

V 

,-68.98877d0 

,3*0.d0 

e 

,89.999998317d0 

,-.5567500297d0 

,1.16119d-5 

$ 

,1.185607d-4 

a 

,53.0066d0 

,-.0566d0 

$ 

,2*0.d0 


j 

,64.5528d0 

,3*0.d0 

8 

,83.31049d0 

,3*0.d0 

U 

,14.920d0 

,3*0.d0 

n 

,41.46635d0 

,3*0.d0 

P 

,66.36420d0 

,3*0.d0 

z 

,52*90.d0 


z 

/ 



c pole vector "i'm not sure ... but it's needed!” 
c 

data ((polev(ij,3),i=l,4) j=0,22) 
s /60.38908d0,-. 10213d0,.00003d0,0.d0, 
m 37.95923d0,-.09577d0,-.00008d0,0.d0, 
v 180.28229d0,.11043d0,.00062d0,0.d0, 
e 89.9999986565d0,-.6402781d0,-3.042075d-4,-5.0837d-6, 
a 43.34526d0,-.09181d0,.0001d0,0.d0, 
j 317.92374d0,.08006d0,-.00019d0,0.d0, 

6 46.0629d0,.01624d0,-.00010d0,0.d0, 
u 6.05795d0,.00182d0,-.00004d0,0.d0, 
n 17.21185d0,.00078d0,-.00002d0,0.d0, 
p 69.13587d0,0.d0,0.d0,0.d0, 

$ 52*0.d0/ 
c 

c *** infvc common 
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data esnprt / O.dO / 
data pine, pinctm / 1.0d+6, O.dO/ 
data pmc, pmctm /-l.dO, O.dO / 
c 

c include 'motblc.inc' 


c common /motblc / 

c motblc is initialized in subroutine nomtab 
data 

$ thrstt(l) /0. / 

c 

c *** paeve common 
c 

data pacav / 3*0.d0/ 


c 

c 

c 


include 'phzvc.inc' 


c common /phzvc / 

data 

$ maxtim /1000. / 

c 

c *** common planec 
c 

data altatm / O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, 

$ O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, 0.d0, 

$ O.dO, O.dO, O.dO, O.dO, O.dO, O.dO, O.dO / 

data gj2 /0.d0,.0001d0,0.d0,.00108263d0,.0019628d0„014733d0 

$ 01664d0,.023d0,0.d0,0.d0,.00021d0,12*0.d0/ 

data re /69600 () . d0 , 2439 . d0 , 6()52 . d0 > 6 3 78 i 4 d 0 | 3397 2d0,71398.d0, 

$ 60000.d0,25400.d0,24300.d0,1500.d0,1738.09d0,12*0.d0/ 

c 

data thedot /2.86533d-6, 1.2400 12d-6^.992369d-7,7.292116d-5 
$ 7.088218d-5,1.758532d-4,1.63785d-4,1.012372d-4, 

$ 9.469017d-5,1.138585d-5,13*0.d0/ 


c 



c common propdc 
c 

data bthrst / I.d0,2*0.d0 / 
data cd / O.dO / 
data d / O.dO / 
datafbody / 2*0.d0/ 
data trad /O.dO/ 
data frac /O.dO/ 
data go /9.80665d-3/ 
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datapdump /O.dO / 
data pitchO / O.dO / 
data pitchl / O.dO / 
data pitch2 / O.dO / 
data prpdat / 20*0.d0 / 
data pscale / l.dO / 
datarect /.OldO/ 
datarollO /O.dO/ 
datarolll /O.dO/ 
data roll2 /O.dO/ 
data scdrag / 3*0.d0/ 
data scmass / l.d6 / 
data scsfa / 3*0.d0/ 
data sfc / 1.024dl7/ 
data spi / l.d6/ 
data step /O.dO/ 
data thnist / l.d6 / 
data thlO /l.dO/ 
data thll /O.dO/ 
data thl2 /O.dO/ 
data wjett / O.dO / 
data wjettm / O.dO / 
data wprop / l.d6 / 
data yawO /O.dO/ 
data yawl / O.dO / 
data yaw2 / O.dO / 
c 

c *** sensit common 
c 

data pertsj / 2*l.d-6 / 
c 

c *** statec common 
c 

data x / 6*0.d0/ 
data xl / 6*0.d0/ 


c 

c 

c *** tgovc common 
c 

data fuxn /6hfuxnl ,6hfuxn2 ,6hfuxn3 ,6hfuxn4 ,6hfuxn5 , 
$ 6hfuxn6 ,6hfiaxn7 ,6hfUxn8 ,6hfuxn9 ,6hfuxnl0/ 

c 

data tgo /6htgo / 
data timemn/6htimemn / 
data timx /6htimx / 
c 

end 
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block data datal 
implicit integer (a-z) 

c 

c datal - defines computational commons 
c (aaiv - zzend) computational values 
c 

c the data statements in this routine are standardized 
c 

C ***:M>******************************************* 

c 

include 'collonc.inc' 
include 'cycnc.inc' 
include 'ephejc.inc' 
indude 'infvlc.inc' 
include 'mnvrjc.inc' 
indude 'phznc.inc' 
indude 'propjc.inc' 
indude 'statnc.inc' 
include 'tgonc.inc' 
include 'zzendc.inc' 

c common /collonc/ 

data nsegph /V 
data nsgphO /0/ 
data nsgpwd /-l/ 
data nsgpwi /-V 

c common /collocc/ 

data cointp / linear'/ 

c common /cycnc / 

data 

1 jevt 10 / 

c *** cyccc common 
c 

data iepoch / ’julian'/ 
data iprop / 'lstep'/ 
c 

c *** ephecc common 
c 

data kernel //home/clarkent/pephem/leapseconds .ker ' , 
$ '/home/darkent/pephem/p_constants.ker7 

data pefile f condensed'/ 
data planet f sun ', 
m 'mercury ', 

v Venus 

e 'earth ', 

m 'mars ', 
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9 
t 

9 
I 

9 

7 

data spk f /home/clarkent/pephem/de 130.bsp’,9*'none V 
c 

c *** infvlc common 
c 

data title / ' V 

data fid / ’ipost file i.d. 00007 
data sfid / 
c 

c *** mnvrjc common 
c 

data ilnch /O/ 
c 

c *** mnvrcc common 
c 

data mantyp P non e7 
c 

c *** phznc common 
c 

data inff / 0/ 
data ipif / 1/ 
c 

c *** propjc common 
c 

data idt / 1 / 
data icoord 111 
data itabsz 122 / 

data npert / 0,0, 0,0, 0,0, 0,0,0, 0,0,0, 0,0, 0,0,0, 0,0,0, 0/ 
data iforce / 0,0, 0,0,0 / 
data iatmos 111 
datajett / 0 / 
c 

c *** propcc common 
c 

data massfi / 'none' / 
c 

c *** statue common 
c 

dataidbody/O/ 


'jupiter ’, 
'saturn ’, 
'uranus ', 
'neptime ’, 
'pluto ’, 
'moon ', 

t t i t t 


J 
s 
u 
n 

P 
m 
@ 

@ 

data ref A>1950'/ 


» 

i r 


t i 
9 

1 f 
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data ipbody / 0,3/ 

data ibody /0000, 0001, 0002, 0003, 0004, 0005, 0006, 0007, 0008, 0009, 
+ 0103,12*0/ 

data jplbdy/10, 1,2, 399, 499, 599, 699, 799, 899, 999, 301, 12*0/ 
datajpbody/-l/ 
c 

c statcc common 
c 

data idfram / 'ecliptic', 'mean2000'/ 
data inputx / 'none'/ 
c 
c 

c include 'tgonc.inc' 

c 

c common /tgonc / 

data 

1 nxevt /0 / 

c 

c include 'zzendc.inc' 



c common /zzendc / 

data 

1 end /0. / 

c 

end 
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block data diet 
implicit real*8 (a-h,o-z) 

c *********************************************************** 

c 

c diet - defines computational commons 
c (aaiv - zzend) dictionary values data 
c 

c the data statements in this routine are standarized. 
c 

c *********************************************************** 

c 

include 'aaivd.inc' 
include 'atmosd.inc' 
include 'cycvd.inc' 
include 'dvivd.inc' 
include 'ephemd.inc' 
include 'infvd.inc' 
include 'motbld.inc' 
include 'pacvd.inc' 
include 'phzvd.inc' 
include 'planed.inc' 
include 'propdd.inc' 
include 'sensid.inc' 
include 'stated.inc' 
include 'tgovd.inc' 
c 

c *********************************************************** 

C 

c *** commons not included in dictionary 
c 

£*♦♦♦♦♦*♦*♦♦*♦+♦***♦**♦♦♦*♦♦*****♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦♦ 

c 

c none, 
c 

c ************************************************ *********** 
c 

c *** computational data dictionary 
c 

c *********************************************************** 

c 

c *** aaivd common 
c 

data vi /6hvi0 ,6hvil ,6hvi2 / 
c 

c *** atmosd common 
c 

data alt /6haltl ,6halt2 ,6halt3 ,6halt4 ,6halt5 , 

+ 6halt6 ,6halt7 ,6halt8 ,6halt9 ,6haltl0 , 

+ 6haltll ,6haltl2 ,6haltl3 ,6haltl4 ,6haltl5 , 
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+ 6haltl6 ,6haltl7 ,6haltl8 ,6haltl9 ,6halt20 , 

+ 6halt21 ,6halt22 ,6halt23 ,6halt24 ,6halt25 , 

+ 6halt26 ,6halt27 ,6halt28 .6halt29 ,6halt30 , 

+ 6halt31 ,6halt32 ,6halt33 • balt34 ,6halt35 , 

+ 6halt36 ,6halt37 ,6halt38 , ialt39 ,6halt40 , 

+ 6halt41 ,6halt42 ,6halt43 ,onalt44 ,6halt45 , 

+ 6halt46 ,6halt47 ,6halt48 ,6halt49 ,6halt50 / 

data pres /6hpresl ,6hpres2 ,6hpres3 ,6hpres4 ,6hpres5 , 

+ 6hpres6 ,6hpres7 ,6hpres8 ,6hpres9 ,6hpresl0, 

+ 6hpresll,6hpresl2,6hpresl3,6hpresl4,6hpre8l5, 

+ 6hpresl6,6hpresl7,6hpresl8,6hpresl9,6hpres20, 

+ 6hpres2 1 ,6hpres22 ,6hpres23 , 6hpres24 ,6hpres25 , 

+ 6hpres26,6hpres27,6hpres28,6hpres29,6hpres30, 

+ 6hpres31,6hpres32,6hpres33,6hpres34,6hpre&35, 

+ 6hpres36,6hpres37,6hpres38 > 6hpres39,6hpres40, 

+ 6hpres41,6hpres42,6hpres43,6hpres44,6hpres45, 

+ 6hpres46,6hpres47 > 6hpres48,6hpres49,6hpres50 / 

datatpt /6htptl ,6htpt2 ,6htpt3 ,6htpt4 ,6htpt5 , 

+ 6htpt6 ,6htpt7 ,6htpt8 ,6htpt9 ,6htptl0 , 

+ 6htptll ,6htptl2 ,6htptl3 ,6htptl4 ,6htptl5 , 

+ 6htptl6 ,6htptl7 ,6htptl8 ,6htptl9 ,6htpt20 , 

+ 6htpt21 ,6htpt22 ,6htpt23 ,6htpt24 ,6htpt25 , 

+ 6htpt26 ,6htpt27 ,6htpt28 ,6htpt29 ,6htpt30 , 

+ 6htpt31 ,6htpt32 ,6htpt33 ,6htpt34 ,6htpt35 , 

+ 6htpt36 ,6htpt37 ,6htpt38 ,6htpt39 ,6htpt40 , 

+ 6htpt41 ,6htpt42 ,6htpt43 ,6htpt44 ,6htpt45 , 

+ 6htpt46 ,6htpM7 ,6htpt48 ,6htpt49 ,6htpt50 / 
data wgt /6hwgtl ,6hwgt2 ,6hwgt3 ,6hwgt4 ,6hwgt5 , 

+ 6hwgt6 ,6hwgt7 ,6hwgt8 ,6hwgt9 ,6hwgtl0 , 

+ 6hwgtll ,6hwgtl2 ,6hwgtl3 ,6hwgtl4 ,6hwgtl5 , 

+ 6hwgtl6 ,6hwgtl7 ,6hwgtl8 ,6hwgtl9 ,6hwgt20 , 

+ 6hwgt21 ,6hwgt22 ,6hwgt23 ,6hwgt24 ,6hwgt25 , 

+ 6hwgt26 ,6hwgt27 ,6hwgt28 ,6hwgt29 ,6hwgt30 , 

+ 6hwgt31 ,6hwgt32 ,6hwgt33 ,6hwgt34 ,6hwgt35 , 

+ 6hwgt36 ,6hwgt37 ,6hwgt38 ,6hwgt39 ,6hwgt40 , 

+ 6hwgt41 ,6hwgt42 ,6hwgt43 ,6hwgt44 ,6hwgt45 , 

+ 6hwgt46 ,6hwgt47 ,6hwgt48 ,6hwgt49 ,6hwgt50 / 

data atms /6hatmsl ,6hatms2 ,6hatms3 ,6hatms4 ,6hatms5 , 
+ 6hatms6 ,6hatms7 ,6hatms8 ,6hatms9 ,6hatmsl0/ 

c 

c *** cycvd common 
c 

data date 

1 /6hyear ,6hmonth ,6hday ,6hhour ,6hminute,6hsecond/ 
data dt /6hdt / 
data dtime /6hdtime / 

data dtimr /6hdtimr 1 ,6hdtimr2 ,6hdtimr3 ,6hdtimr4 , 6hdtimr5 , 
+ 6hdtimr€,6hdtimr7,6hdtiinr8,6hdtinir9,6hdtiml0, 

+ 6hdtim 1 1 ,6hdtiml2, 6hdtim 1 3 ,6hdtim 14 ,6hdtim 15 , 
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+ 6hdtiml6,6hdtiml7,6hdtiml8,6hdtiml9,6hdtim20/ 

data dtO /6hdt0 / 
data juldat/6hjuldat/ 
data tdurp /6htdurp / 
data time /6htime / 

data timrb /6htimrbl,6htimrb2,6htimrb3,6htimrb4,6htimrb5 ) 
+ 6htimrb6,6htimrb7,6htimrb8,6htimrb9,6htimbl0, 

+ 6htimbll,6htimbl2 > 6htimbl3 > 6htimbl4,6htiinbl5, 

+ 6htimbl6,6htimbl7,6htimbl8,6htiinbl9,6htimb20/ 

data timrf /6htimrfl,6htimrf2 ,6htimrf3 ,6htimrf4,6htimrf5 , 

+ 6htimrfB,6htimrf7,6htimrfB,6htiinrf9,6htiinfl0, 

+ 6htimfl I,6htimn2,6htimfl3,6htimfl4,6htimfl5, 

+ 6htimfl6,6htimfi7.6htimfl8.6htimfi9,fih timf9 n/ 

data tref /6htref / 
data trefp /6htrefp / 
c 

c *** dvivd common 


data dvx /6hdvx ,6hdvy ,6hdvz / 
data dvmag /6hdvmag / 
data dvsum /6hdvsum / 


c 

c *** ephemd common 
c 


c 


c 


c 


data gmu /6hgmu0 ,6hgmul ,6hgmu2 ,6hgmu3 ,6hgmu4 
+ 6hgmu5 ,6hgmu6 ,6hgmu7 ,6hgmu8 ,6hgmu9 , 

+ 6hgmulO ,6hgmull ,6hgmul2 ,6hgmul3 ,6hgmul4 , 

+ 6hgmul5 ,6hgmul6 ,6hgmul7 ,6hgmul8 ,6hgmul9 , 

+ 6hgmu20 ,6hgmu21 ,6hgmu22 / 


data rsoi /6hrsoi0 ,6hrsoil ,6hrsoi2 ,6hrsoi3 ,6hrsoi4 , 
+ 6hrsoi5 ,6hrsoi6 ,6hrsoi7 ,6hrsoi8 ,6hrsoi9 , 

+ 6hrsoi 1 0 ,6hrsoi 1 l,6hrsoil2,6hrsoil3,6hrsoil4, 

+ 6hrsoil5,6hrsoil6,6hrsoil7,6hrsoil8,6hrsoil9, 

+ 6hrsoi20,6hrsoi2 l,6hrsoi22/ 


data tsoi /6htsoi0 ,6htsoil ,6htsoi2 ,6htsoi3 ,6htsoi4 , 
+ 6htsoi5 ,6htsoi6 ,6htsoi7 ,6htsoi8 ,6htsoi9 , 

+ 6htsoi 10 ,6htsoi 1 1 ,6htsoi 12 ,6htsoi 13,6htsoi 14, 

+ 6htsoil5,6htsoil6,6htsoil7,6htsoil8,6htsoil9, 

+ 6htsoi20 ,6htsoi2 1 ,6htsoi22/ 


data (elem(i),i=l,50) 

^ /6heleml ,6helem2 ,6helem3 ,6helem4 ,6helem5 , 

f 6helem6 ,6helem7 ,6helem8 ,6helem9 ,6heleml0, 

h 6helem 1 1 ,6helem 1 2,6helem 13 ,6heleml4 ,6helem 15 , 

h 6heleml6,6heleml7,6heleml8,6heleml9,6helem20, 

h 6helem21,6helem22,6helem23,6helem24,6helem25, 

h 6helem26,6helem27,6helem28,6helem29,6helem30, 
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6helem3 1 ,6helem32,6helem33,6helem34,6helem35 , 
6helem36,6helem37,6helem38,6helem39,6helem40, 
6helem4 1 ,6helem42,6helem43,6helem44,6helem45 , 
6helem46,6helem47,6helem48,6helem49,6helem50/ 


data (elem(i)4=51,100) 

+ /6helem5 1 ,6helem52,6helem53 ,6helem54,6helem55 , 

+ 6helem56 ,6helem5 7 ,6helem58 ,6helem59,6helem60 , 

+ 6helem61,6helem62,6helem63,6helem64,6helem65, 

+ 6helem66,6helem67,6helem68,6helem69,6helem70, 

+ 6helem71,6helem72,6helem73,6helem74,6helem75, 

+ 6helem76,6helem77,6helem78,6helem79,6helem80, 

+ 6helem81,6helem82,6helem83,6helem84,6helem85, 

+ 6helem8P.6helem87,6helem88,6helem89,6helem90, 

+ 6helem9: , 6helem92 ,6helem93 ,6helem94,6helem95 , 

+ 6helem96,6helem97,6helem98,6helem99,6helel00/ 

data (elem(i),i=10 1,150) 

+ /6helel01,6helel02,6helel03,6helel04,6helel05, 

+ 6helel06,6helel07,6helel08,6helel09,6helell0, 

+ 6helelll,6helell2,6helell3,6helell4,6helell5, 

+ 6helell6,6helell7,6helell8,6helell9,6helel20, 

+ 6helel21,6helel22,6helel23,6helel24,6helel25, 

+ 6helel26,6helel27,6helel28,6helel29,6helel30, 

+ 6helel31,6helel32,6helel33,6helel34,6helel35, 

+ 6helel36,6helel37,6helel38,6helel39,6helel40, 

+ 6helel41,6helel42,6helel43,6helel44,6helel45, 

+ 6hele 146,6hele 14 7 ,6hele 148 ,6hele 149,6helel50/ 


data (elem(i),i=151,200) 

+ /6helel51,6helel52,6helel53,6helel54,6helel55, 

+ 6helel56,6helel57,6helel58,6helel59,6hdel60, 

+ 6helel61,6helel62,6helel63,6helel64,6helel65, 

+ 6helel66,6helel67,6helel68,6helel69,6helel70, 

+ 6helel71,6helel72,6helel73,6helel74,6helel75, 

+ 6helel76,6helel77,6helel78,6helel79,6helel80, 

+ 6helel81,6helel82,6helel83,6helel84,6helel85, 

+ 6helel86,6helel87,6helel88,6helel89,6helel90, 

+ 6helel91,6helel92,6helel93,6helel94,6helel95, 

+ 6helel96,6helel97,6helel98,6helel99,6hele200/ 

data (elem(i),i=201,250) 

+ /6hele201,6hele202,6hele203,6hele204,6hele205, 

+ 6hele206,6hele207,6hele208,6hele209,6hele210, 

+ 6hele21 1,6hele212,6hele213,6hele214,6hele215, 

+ 6hele216,6hele217,6hele218,6hele219,6hele220, 

+ 6hele22 l,6hele222 ,6hele223 ,6hele224 ,6hele225 , 

+ 6hele226,6hele227,6hele228,6hele229,6hele230, 

+ 6hele231,6hele232,6hele233,6hele234,6hele235, 

+ 6hele236,6hele237,6hele238,6hele239,6hele240, 

+ 6hele241,6hele242,6hele243,6hele244,6hele245, 

+ 6hele246,6hele247,6hele248,6hele249,6hele250/ 

data (elem(i),i=25 1,300) 
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+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


/6hde251,6hde252,6hde253,6hde254,6hele255, 
6hde256,6hde257,6hde258,6hde259,6hde260, 
6hde261,6hde262,6hde263,6hde264,6hele265, 
6hde266,6hde267,6hde268,6hde269,6hde270, 
6hde27 1,6hde272,6hde273,6hele274,6hde275, 
6hde276,6hde277,6hde278,6hde279,6hde280, 

6hele281,6hele282,6hele283,6hele284,6hele285, 

6hele286,6hele287,6hele288,6hele289,6hele290, 

6hele291,6hele292,6hele293,6hele294,6hele295, 

6hele296,6hele297,6hele298,6hele299,6hele300/ 


data (elem(i),i = 301350) 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


/6hele301,6hele302,6hele303,6hele304,6hele305, 
6hele306,6hele307,6hele308,6hele309,6hele310, 
6hde311,6hde312,6hde313,6hde314,6hele315, 
6hele316,6hele317,6hele318,6hele319,6hde320, 
6hele321,6hele322,6hele323,6hele324,6hele325, 
6hele326,6hele327,6hele328,6hele329,6hele330, 
6hde331,6hde332,6hde333,6hde334,6hele335, 
6hde336,6hde337,6hde338,6hde339,6hele340, 
6hele34 1 ,6hele34 2 ,6hele34 3 ,6hele344 ,6hele345, 
6hde346,6hde347,6hde348,6hde349,6hde350/ 


data (elem(i),i=35 1,400) 

/6hele35 1,6hele352,6hele353 ,6hele354,6hele355 , 
6hele35 6 ,6hele35 7 ,6hele35 8 ,6hele359,6hele360 , 
6hele36 1 ,6hele362 ,6hele363 ,6hele364 ,6hele365 , 

6hde366,6hde367,6hde368,6hde369,6hde370, 
6hde371,6hde372,6hde373,6hde374,6hde375, 
6hde376,6hde377,6hde378,6hde379,6hde380, 
6hde381,6hde382,6hde383,6hde384,6hde385, 

6hele386 ,6hele38 7 ,6hele38 8 ,6hele389 ,6hele390 , 
T 6hde391,6hde392,6hde393,6hde394,6hde395, 

+ 6hele396,6hele397,6hele398,6hele399,6hele400/ 

dat a (dem(i),i=401,450) 

/6hele40 1 ,6hele402 ,6hele403 ,6hele404,6hele405 , 
6hele40 6 ,6hele40 7 , 6hele40 8 ,6hele409,6hele4 1 0 , 
6hele4 1 1 ,6hele4 1 2 ,6hele4 1 3 ,6hele4 14,6hde4 1 5 , 
6hde416,6hde417,6hde418,6hde419,6hde420, 
6hde421,6hde422,6hde423,6hde424,6hde425, 
6hde426,6hde427,6hde428,6hde429,6hde430, 
6hele431,6hele432,6hde433,6hde434,6hde435, 
6hele436,6hele437,6hde438,6hde439,6hde440, 
6hde441,6hde442,6hde443,6hde444,6hele445, 
6hde446,6hde447,6hde448,6hde449,6hde450/ 

data (elem(i),i=45 1,500) 

+ /6hde451,6hde452,6hele453,6hde454,6hele455, 

+ 6hde456,6hde457,6hde458,6hde459,6hde460, 

+ 6hele461,6hde462,6hele463,6hde464,6hele465, 

+ 6hele466,6hde467,6hde468,6hele469,6hde470, 

+ 6hele471,6hele472,6hde473,6hde474,6hde475, 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 
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+ 

+ 

+ 

+ 

+ 

+ 

+ 


+ 6hele47 6,6hele477,6hele478,6hele479,6hele480 

+ 6hele48 I,6hde482,6hde483,6hde484,6hde485 

+ 6hele486,6hele487,6hele488,6hele489,6hele490, 

+ 6hele491,6hele492,6hele493,6hele494,6hele495, 

+ 6hele496,6hele497 ,6hde498,6hde499,6hde500/ 

data (elem(i),i=50 1,550) 

+ /6hele501,6hele502,6hele503,6hele504,6hele505 

+ 6hele506,6hele507,6hele508,6hele509,6hele510 

6hele5 1 1 ,6hde5 12,6hde5 13,6hde5 14,6hde515 ’ 
6hde516,6hde517,6hde518,6hde519,6hde520 
6hde521,6hde522,6hde523,6hde524,6hde525 
6hele526,6hde527,6hde528,6hde529,6hde53o’ 
6hele531,6hele532,6hele533,6hde534,6hde535’ 
6hele536,6hde537,6hde538,6hde539,6hde54o’ 
6hde541,6hde542,6hde543,6hde544,6hde545’ 

, . , 6hele546 6hde547,6hde548,6hde549,6hele550/ 

data (elem(i), i=551, 552) 

+ /6hele551,6hele552/ 

data (pole(i)4=l,50) 

+ /Ihpojel ,6hpde2 ,6hpde3 ,6hpole4 ,6hpole5 , 

6hjx)l e 6 6hp°le7 ,6hpde8 ,6hpde9 ,6hpolel0, 

fk P °i e ^l , fu PO J el2 ’ 6hpolel3,6h P olel4 » 6h P°*el5, 

fi? PO i e oi , fu P °! e22,6hpole23 » 6h P° le24 * 6h Pole25, 
6h p ol e 26,6h po l e27 , 6hpole28 > 6h p° ie29 

I^ P °? e fi’f^ PO ? e 22' 6 f lpole33 ' 6h P ole34 > 6h Pole35, 

fiK PO } e f?,cu PO i e37 ,^ h Po le 3 8 , 6h P ole 39,6hpole40, 

. fiu PO i 1 !^i’fu POe42,6hpole43 '6hpole44,6hpole45 

data(p 0 fXluX e4,6hPOl648 ' 6hp0le49 ' 6hp ° ,eW 

+ /6hpol e 51,6hpol e 52,6hpole53 ) 6hpole54,6hpole55, 

6h poe ^,6h poe67 , 6h p ole 6 8 , 6 h P o l e6 9 >6h ^ 1 e70 
«K P °| e ^’?u PO ! e72,6h ^ le73 > 6h P° le74 .6hpole75 

«J? P °} e I^fl lp ^ e77 ’ 6hpole78 « 6h P ole79 > 6hpole80 
data(pd*i«S 

+ / ^l 1 ^'f! lpO , 1 , 1 !? 2 ’ 6hpoll03 ’ 6h P oll04 .6hpoll05, 

6h po lM,6h p ol 0 7, 6 hpon°8, 6h p o 1109,6hpolllO, 

6h^?}fi «^}, 1 ?'fu r “! 113 ' 6hpollH ' 6h P“ 111 5 
f k°i ioi^L P °!!F' 6 ! lp0 1 18 ' 6h P°l 1 19,6hpol 120, 

6hpoJ121,6hpoll22,6hpoll23,6hpo]124,6hDoll25 

«^°???^ po ! l27 ' 6h ^ i2 «’ 6h ^ i29 . 6h Mi35: 

6 hpoll 31 , 6 hpoll 32 , 6 hpoll 33 , 6 hpdl 34 , 6 hi>oll 35 , 


+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 


+ 

+ 

+ 

♦ 

+ 

+ 

♦ 

+ 


+ 

+ 

+ 


+ 

+ 
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/6hpoll51,6hpoll52,6hpoll53,6hpoll54,6hpoll55, 
6hpol 156,6hpol 157 ,6hpol 158,6hpoll59,6hpoll60 , 

6 hpoll 61 , 6 hpoll 62 , 6 hpoll 63 , 6 hpoll 64 , 6 hpoll 65 , 

6hpol 166,6hpol 167 ,6hpol 168,6hpoll69,6hpoll70 , 
6hpoll71,6hpoll72,6hpoll73,6hpoll74,6hpoll75, 
6hpol 176,6hpol 177 ,6hpol 178 ,6hpoll79,6hpoll80, 
6hpoll81,6hpoll82,6hpoll83,6hpoll84,6hpoll85, 
6hpoll86,6hpoll87,6hpoll88,6hpoll89,6hpoll90, 
6hpoll91,6hpoll92,6hpoll93,6hpoll94,6hpoll95, 
6hpol 196,6hpol 197 ,6hpol 198,6hpoll99,6hpol20Q/ 


+ 6hpol 136 ,6hpol 137 ,6hpol 138 ,6hpoll39,6hpol 140 , 

+ 6hpoll41,6hpoll42,6hpoll43,6hpoll44,6hpoll45, 

+ 6hpol 146,6hpol 147 ,6hpol 148 ,6hpol 149,6hpoll50/ 

data (pole(i),i=151,200) 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

+ 

data (pole(i)4=201, 250) 

+ /6hpol20 1 ,6hpol202 ,6hpol203 ,6hpol204,6hpol205, 

+ 6hpol206,6hpol207,6hpol208,6hpol209,6hpol210, 

+ 6hpol211,6hpol212,6hpol213,6hpol214,6hpol215, 

+ 6hpol2 1 6 , 6hpol2 1 7 , 6hpol2 18 , 6hpol2 19 ,6hpol220 , 

+ 6hpol22 1 ,6hpol222 ,6hpol223,6hpol224,6hpol225, 

+ 6hpol226,6hpol227,6hpol228,6hpol229,6hpol230, 

+ 6hpol231,6hpol232,6hpol233,6hpol234,6hpol235, 

+ 6hpol236,6hpol237,6hpol238,6hpol239,6hpol240, 

+ 6hpol241,6hpol242,6hpol243,6hpol244,6hpol245, 

+ 6hpol246,6hpol247,6hpol248,6hpol249,6hpol250/ 

data (pole(i)4=25 1,276) 

+ /6hpol251,6hpol252,6hpol253,6hpol254,6hpol255, 

+ 6hpol256,6hpol257 ,6hpol258 ,6hpol259,6hpol260, 

+ 6hpol261,6hpol262,6hpol263,6hpol264,6hpol265, 

+ 6hpol266,6hpol267,6hpol268,6hpol269,6hpol270, 

+ 6hpol271,6hpol272,6hpol273,6hpol274,6hpol275, 

+ 6hpol276/ 


c *** infvd common 


data esnprt / 6hesnprt/ 
data pine, pinctm / 6hpinc ,6hpinctm/ 
data pmc, pmctm / 6hpmc ,6hpmctm/ 
c 

c *** motbld common 


data thrstt / 6hthrstt/ 
data cdt / 6hcdt / 
data dt / 6hclt / 
data fmasit / 6hfinasit/ 
data finasst / 6hfinasst/ 
data smasst / 6hsmasst/ 


c 

c 
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c include 'pacvd.inc' 


c common pacvd 
data 

1 pacav /6hpacavx,6hpacavy,6hpacavz/ 
c 

c *** phzvd common 
c 

data maxtim/6hmaxtim/ 
data tcpmax/6htcpmax/ 
c 


c *** planed common 
c 

data altatm/6haltat0,6haltatl,6haltat2,6haltat3,6haltat4, 
+ 6haltat5,6haltat6,6haltat7,6haltat8,6haltat9, 

+ 6haltal0,6haltal I,6haltal2,6haltal3,6haltal4, 

+ 6haltal5,6haltal6,6haltal7,6haltal8,6haltal9, 

+ 6halta20,6halta21,6halta22/ 

data gj2 /6hgj20 ,6hgj21 ,6hgj22 ,6hgj23 ,6hgj24 , 

+ 6hgj25 ,6hgj26 ,6hgj27 ,6hgj28 ,6hgj29 , 

+ 6hgj210 ,6hgj211 ,6hgj212 ,6hgj213 ,6hgj214 , 

+ 6hgj215 ,6hgj216 ,6hgj217 ,6hgj218 ,6hgj219 , 

+ 6hgj220 ,6hgj221 ,6hgj222 / 

data re /6hre0 ,6hrel ,6hre2 ,6hre3 ,6hre4 , 

+ 6hre5 ,6hre6 ,6hre7 ,6hre8 ,6hre9 , 

+ 6hrel0 ,6hrell ,6hrel2 ,6hrel3 ,6hrel4 , 

+ 6hrel5 ,6hrel6 ,6hrel7 ,6hrel8 ,6hrel9 , 

+ 6hre20 ,6hre21 ,6hre22 / 

data thedot/6hthdt0 ,6hthdtl ,6hthdt2 ,6hthdt3 ,6hthdt4 
+ 6hthdt5 ,6hthdt6 ,6hthdt7 ,6hthdt8 ,6hthdt9 , 

+ 6hthdt i 0 , 6 hthdt i 1 , 6 hthdt 3.2 , 6hthdt 13 ,6hthdtl4 , 

+ 6hthdtl5,6hthdtl6,6hthdtl7,6hthdtl8,6hthdtl9, 

+ 6hthdt20 ,6hthdt2 1 ,6hthdt22/ 

c 

c *** propdd common 
c 


data bthrst /6hbthrsl > 6hbthrs2,6hbthrs3/ 

data cd /6hcd / 

data cl /6hd / 

data fbody /6hfbodyl,6hfbody2/ 

data frad /6hfraci / 

data frac /6hfrac / 

data go /6hgo / 

data pdump /6hpdump/ 

data pitchO /6hpitch0/ 

data pitchl /6hpitchl/ 

data pitch2 /6hpitch2/ 

data pt /6hpitch / 

data prpdat /6hprpdl ,6hprpd2 ,6hprpd3 ,6hprpd4 ,6hprpd5 , 
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h 6hprpd6 ,6hprpd7 ,6hprpd8 ,6hprpd9 ,6hpipdlO, 

h 6hprpdll,6hprpdl2,6hprpdl3,6hprpdl4,6hprpdl5, 

t- 6hprp d 1 6 , 6hprpd 1 7 , 6hprpd 1 8 , 6hprpd 19 ,6hprpd2 0/ 

data pscale /6hpscale/ 

datarect /6hrect / 

data rl /6hroll / 

data rollO /6hroll0 / 

data rolll /6hrolll / 

data roll2 /6hroll2 / 

data rose /6hrosc / 

da ta sedrag /6hscdrgl,6hscdrg2,6hscdrg3/ 

data semass /6hscmass/ 

data sesfa /6hscsfal,6hscsfa2,6hscsfa3/ 

data sfc /6hsfc / 

data spi /6hspi / 

data step /6hstep / 

data thrust /6hthrust/ 

data thlO /6hthl0 / 

data thll /6hthl 1 / 

data thl2 /6hthl2 / 


datatprev /6htprev/ 
data wjett /6hwjett / 
data wjettm /6hwjettm/ 
data wprop /6hwprop / 

data xosc /6hxoscl ,6hxosc2 ,6hxosc3 ,6hxosc4 ,6hxosc5 , 
+ 6hxosc6 / 

data xtrue /6hxtrul ,6hxtru2 ,6hxtru3 ,6hxtru4 ,6hxtru5 , 
+ 6hxtru6 ,6hxtru7 ,6hxtru8 ,6hxtru9 ,6hxtrul0, 

+ 6hxtrull,6hxtrul2,6hxtrul3,6hxtrul4,6hxtrul5, 

+ 6hxtrul6,6hxtrul7,6hxtrul8,6hxtrul9,6hxtru20, 

+ 6hxtru2 1 ,6hxtru22,6hxtru23,6hxtru24,6hxtru25 , 

+ 6hxtru26,6hxtru27,6hxtru28,6hxtru29,6hxtru30 , 

+ 6hxtru31,6hxtru32,6hxtru33,6hxtru34,6hxtru35, 

+ 6hxtru36 ,6hxtru37,6hxtru38 ,6hxtru39,6hxtru40 , 


f 6hxtru41,6hxtru42/ 

data yawO /6hyaw0 / 
data yawl /6hyawl / 
data yaw2 /6hyaw2 / 
data yw /6hyaw / 


c 

c *** sensid common 


data delta /6hdltll ,6hdltl2 ,6hdltl3 ,6hdltl4 ,6hdltl5 , 
+ 6hdltl6 ,6hdltl7 ,6hdltl8 ,6hdltl9 ,6hdltll0, 

+ 6hdlt21 ,6hdlt22 ,6hdlt23 ,6hdlt24 ,6hdlt25 , 

+ 6hdlt26 ,6hdlt27 ,6hdlt28 ,6hdlt29 ,6hdlt210, 

+ 6hdlt31 ,6hdlt32 ,6hdlt33 ,6hdlt34 ,6hdlt35 , 

+ 6hdlt36 ,6hdlt37 ,6hdlt38 ,6hdlt39 ,6hdlt310/ 
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data epst /6hepst / 

data exvel /6hexvll ,6hexvl2 ,6hexvl3 ,6hexvl4 ,6hexvl5 , 

+ 6hexvl6 ,6hexvl7 ,6hexvl8 ,6hexvl9 ,6hexvll0/ 

data wmass 

1 /6hwmassl,6hwmass2,6hwmass3,6hwmass4,6hwmass5 

2 ,6hwma8s6,6hwmass7,6hwma8s8,6hwmas89,6hwma8l0 

2/ 


c 

c *** stated common 


c 


c 


c 


c 


data orbp / 

1 6hsma , 6heccen , 6hinc f 6hanlong, 6hargp , 

2 6hmeaan , 6htruan , 6hrperi , 6hvperi , 6htfp , 

3 6hrapoap, 6hperiod, 6hradius, 6hspeed , 6hbtheti, 

4 6hvinfiri, 6hvinfyi, ehvinfzi, 6hbdti , 6hbdri , 

5 6hhypta , 6hc3i , 6hrai , 6hded ,6hbtheto, 

6 6hvinfxo, 6hvinfyo, 6hvinfzo, 6hbdto , 6hbdro , 

7 6hhypta , 6hc3o , 6hrao , 6hdeco , 6hfpa , 

6 6haltdt , 6hlongp , 6hlat , 6hlong , 6hb / 

data soi /6hsoi / 
data ti /6hti / 
data tf /6htf / 


datax /5hx ,5hy ,5hz ,5hvx ,5hvy ,5hvz / 
data xeq /5hxeq ,5hyeq ,5hzeq ,5hvxeq ,5hvyeq ,5hvzeq / 
data xo /5hxo ,5hyo ,5hzo ,5hvxo ,5hvyo ,5hvzo / 

j *5hyl ,5hzl ,5hvxl ,5hvyl ,5hvzl / 

data x2 /5hx2 ,5hy2 ,5hz2 ,5hvx2 ,5hvy2 ,5hvz2 / 

data x2eq /5hx2eq ,5hy2eq ,5hz2eq ,5hvx2eq,5hvy2eq,5hvz2eq/ 


c *** tgovd common 
c 


c 


c 


data ftrm /6hfiixnl ,6hfuxn2 ,6hfuxn3 ,6hfuxn4 ,6hfuxn5 
+ 6hfuxn6 ,6hfuxn7 ,6hfuxn8 ,6hfuxn9 ,6hfuxnl0/ 


» 


data tgo /6htgo / 
data timemn/6htimemn / 
data timx /6htimx / 


end 
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block data dictl 
implicit integer (a-z) 

£******************************************##**** 

c 

c dictl • defines computational commons 
c (aaiv - zzend) dictionary values 
c 

c the data statements in this routine are standardized 
c 

c ***** ******************************************* 
c 

include 'collond.inc' 
include ’cycnd.inc' 
include 'ephejd.inc' 
include 'infvld.inc' 
include mnvrjd.inc' 
include 'phznd.inc' 
include 'propjd.inc' 
include 'statnd.inc' 
include 'tgond.inc' 
include 'zzendd.inc' 
c 

c include 'collond.inc' 



c common /collond / 

data nsegph /0/ 
data cointp / linear'/ 
c 



c include ’cycnd.inc' 

c common /cycnd / 

data 

1 jevt /0 / 

c common /cyccd / 

data 

1 iepoch /julian ' / 

c 

c include 'ephecd.inc' 



c common /ephecd / 

data 

1 planet(0) /sun ' / 

c 



c include 'infvld.inc' 
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c common /infvld / 

data 

1 title r / 

c 

c include 'mnvrjd.inc' 

c common /mnrvjd / 

data 

1 ilnch /O / 

c common /mnvrcd / 

data 

1 mantyp /'none' / 

c 

c include 'phznd.inc' 

c 

c common /phznd / 

data 

1 inff /O / 

c 

c include 'propjd.inc' 

c 

c common /propjd / 

data 

1 idt /I / 

data massfi / 'none' / 
c 

c 

c include 'statnd.inc' 

c common /statnd / 

Hfltn 

1 idbody /0 / 

c common /stated / 

1 inputx /’none’ / 

c 

c 

c include 'tgond.inc' 

c common /tgond / 

data 

1 nxevt /O / 

c 

c include 'zzendd.inc' 

c 

c common /zzendd / 
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1 end /O. / 

c 

end 
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program master 
implicit real*8 (a-h,o-z) 
c 

C*** master - main program of the interplanetary post simulation 
c master - main program control 

c 

c external files and unit numbers 
c 2 dictionary names 

c 3 pscan.sav 

c 18 profil 

c 20 computational iv-to-end 

c 31 npinput opened in nlprg 

c 32 npost3d.rst 

c 33 npost3d.out 

c 38 ipo8t.smo, for pscan 
c 

include 'ephejc.inc' 
c 

cc include ’infvc.inc' 
c 

include 'inpvc.inc' 
c 

include 'pagerc.inc' 

include ’parmsc.inc' 

include 'prOfc.inc' 

character*256 arg 

data arg / V / 
c 

c routine added to handle Sun date and time 
preset = 1 
call fdate pdate) 
cc call time (ptime) 

c get profile name from command line 
icnt s iargc() 

if(icnt.ne.O) then 

call getarg ( 1, arg ) 
pmam = arg 

else 

pmam = 'profil' 
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endif 

c... nomtab: set table length into common, set up table nodes 
call nomtab 

c open a scratch file to save the computational iv-to-end 

c 

call savdat 
c 

c open the npsol summary output file for optimization 

c |# f 

open ( unit=32, status =' unknown', file='ipost.rst , 

$ form='formatted’ , access='sequential’) 

c 

c open the npsol summary output file for optimization 

c 

open ( unit=33, status='unknown', file=’ipost.out' , 

$ form='formatted\ access=’sequential’) 

10 continue 
rewind 31 

call readat 
if ( inpef.ge.O ) then 
call tspxm 
endif 

if ( inpef.eq.O ) go to 10 
c 

if (iephem.eq.l) close (unit=7) 
close ( unit=18, status='keep' ) 
dose ( unit=32, status='keep' ) 
dose ( unit=33, status='keep' ) 
c 

stop 

end 
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9.0 INPUT/OTTTPUT TNTRRE’ftrfift 


POST PRE-PROCESSING: 

Initial conditions, Planetary Ephemeris 



IPOST POST-PROCESSING: 
Interpolation, Extrapolation, 
Transformation, Reformatting 

J 

PLOTS, DISPLAYS, STORAGE 
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